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The refractive index of novel organosilica (nano/micro) material is determined using two methods. 
The first method is based on analysis of optical extinction efficiency of organosilica beads versus 
wavelength, which is obtained by a standard laboratory spectrometer. The second method relies 
on the measurable trapping potential of these beads in the focused light beam (laser tweezers). 
Polystyrene beads were used to test these methods, and the determined dispersion curves of refrac- 
tive index values have been found accurate. The refractive index of organosilica beads has been 
determined to range from 1.60-1.51 over the wavelength range of 300-1100 nm. 



INTRODUCTION 



Properties of (nano/micro)materials are shaped up 
by the bulk and structural properties of their build- 
ing blocks, with an example being refractive index 
that determines the material's optical characteristics. 
These optical properties need to be precisely controlled 
in various applications, such as sizing of organosilica 
(nano/micro)beads in the flow cytometry biological as- 
saying. By manipulating the composition of materi- 
als on the nanoscale, new structural properties have 
been produced. For reliable production of desirable 
(nano/micro)materials both bulk and structural proper- 
ties need to be controlled. Since bulk properties of the 



material can vary as a result of nanotcchnological synthe- 
sis, it is important to monitor them in situ. In particular, 
silica-derived novel material, termed organosilica, have 
been synthesized in the form of (nano/micro)spheres, 
with a range of useful qualities that have been re- 
ported 0, [if. Organosilica surface activation have ap- 
peared to be straightforward to bind oligonucleotide 
and peptide sequences for optical encoding and high- 
throughput screening via flow cytometry 0, 0] ■ The opti- 
cal encoding method, however, demands tight margins on 
organosilica microspheres size and refractive index to en- 
sure reliable readout and correct interpretation of their 
scattering intensities in the context of flow cytometry. 
However, it has been noticed that the refractive index of 
organosilica microspheres is highly dependent on synthc- 



2 



sis route calling for a reliable refractometry method 
amenable for in situ implementation. Since the organosil- 
ica synthesis yields spherical particles ranging from 50 
nm to 6 /an, a number of existing refractometry meth- 
ods applicable for bulk materials are no longer suitable 
for refractometry of this material. Among the suitable 
methods, an index matching method has been reported 
in a number of particle refractometry applications [H, 0] • 
In brief, particles in the solution do not scatter, i.e. ren- 
dering the sample clear, if refractive indices of the par- 
ticles and solution are matched. Therefore, the solution 
refractive index read out at the sample maximum trans- 
mittivity provides a measure of the refractive index of the 
particle ensemble. Despite its conceptual simplicity, the 
method is tedious, requiring preparation of the new so- 
lution of predetermined refractive index per each experi- 
mental data point. The availability of suitable solutions 
with high refractive indices is usually limited. 

The majority of existing refractometry methods for 
particles rely on the acquisition of optical scatter data, 
followed by their fitting to the theoretical model from 
which one can infer the particle refractive index, or its 
size, or both. For example, the determination of re- 
fractive index of an ensemble of microspheres have been 
demonstrated using holographic method of recording of 
the optical scatter [5|. The intensity of optical scatter 
has been parameterized in the Rayleigh scattering model 
to perform sizing of nanodiamonds [6|. Differential mi- 
croscopy measurements of the optical scatter have also 
permitted detection of the minute cell organelle swelling, 
i.e. change in the scatterer's size, as a result of induced 
apoptosis 0- 

Scattering spectroscopy provides another approach in 
the determination of particle refractive index in which 
characteristic particle optical extinction efficiency versus 
wavelength is determined d i, M, EI El El- Usu- 
ally, the acquired spectrum of monodisperse microsphere 
samples is fitted to a theoretical model with the refrac- 
tive index and microsphere size as parameters. Using 
this method, the refractive index for a specified wave- 
length range, i.e. refractive index dispersion, can be de- 
termined. The application of the scattering spectroscopy 
technique to refractometr y o f cells, bacteria, and spores 
has been also reported 0,13, EE EH ■ 

In this paper, we demonstrate the application of two 
refractometry methods to determine the refractive in- 
dex of novel organosilica microspheres. Firstly, we re- 
port application of the refined scattering spectroscopy 
method, where we use a commercial laboratory Lamda40 
spectrophotometer (Perkin Elmer) to measure the re- 
fractive index and dispersion of organosilica, following 
well-characterized polystyrene microsphere refractome- 
try. Secondly, we demonstrate application of the recently 
reported method for measurement of the refractive index 
for single microparticles based on optical tweezers 17 1. 
The two methods are compared, and a very good agree- 



ment is achieved. 



EXPERIMENTAL METHOD 
Materials 

Organosilica spheres were produced using a two-step 
process Q. This process involved the hydrolysis and 
subsequent condensation of thiol-based silica precursor. 
Aqueous solutions of 5.94 /jm organosilica microspheres, 
with the addition of sodium dodecyl sulphate for the pre- 
vention of agglomeration, were prepared for spectroscopy 
measurements. Polystyrene microspheres (ProSciTech, 
5.0% w/v in water) of mean diameters of 2.01, 2.09 and 
5.26 fim, with a coefficient of variance (CV, defined as 
the standard deviation divided by the mean), of less 
than 2.4%, were diluted with MilliQ water. For the laser 
tweezers refractometry, 2.09 /im polystyrene and 5.16 /im 
organosilica microspheres were suspended in MilliQ wa- 
ter. Samples were prepared by placing a drop of the sus- 
pension between the microscope slide and cover slip and 
sealing the chamber with high viscosity silicon grease. 



Methods 

Scanning Electron Microscopy (SEM) 

In order to characterize the size distribution of the syn- 
thesized organosilica microspheres, SEM was used. Sam- 
ples were platinum coated and imaged using SEM (JEOL 
JSM 6400F) with accelerating voltage of 5-10 kV. The 
size distribution of the microspheres was obtained by pro- 
cessing the acquired images followed by statistical anal- 
ysis. 



Scattering Spectroscopy 

In scattering spectroscopy, spherical particle extinc- 
tion efficiency verses illumination wavelength (hereafter 
referred to as the scattering spectrum) was recorded. 
Recorded data were fitted to the theoretical model, where 
fitting parameters represented refractive index and size of 
the particle. To determine the best fit, the value, defined 
as the sum of the difference squared (I exp — Itheor) 2 , 
was minimized. Although this method is simple and ef- 
fective, the final best fit can be significantly unalike. This 
is due to the recording configuration, which should be a 
simple plane wave incident on a single spherical particle 
that cast a shadow on a small-size photodetector situated 
at a distance. The attenuation of the incident optical flux 
along the axis of incidence would be entirely due to light 
scattering by the particle, or in other words, the particle 
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terial. Scans were performed at a rate of 480 nm min -1 
with a spectral resolution of 2 nm. Recorded data were 
digitized and transferred to a personal computer for fur- 
ther processing, archiving, and display. The recorded 
data were fitted to the theoretical scattering spectra, as 
detailed in the modelling method in Section . 



Laser Tweezers Refractometry 



FIG. 1: A lamp produces a wave with divergence of 0.4°. This 
is incident on a quartz cuvette which is filled with an aqueous 
suspension of microspheres. A silicon photodetector measures 
transmission intensity for wavelengths 300-1100 nm. 



Fiber Laser 



scattered light away from the photodetector. Under typ- 
ical experimental conditions, however, illumination from 
an extended light source could not be regarded as the 
plane wave illumination. In addition, instead of a single 
particle, one has to account for an ensemble of particles, 
and a finite photodetector size. 

In this paper, we employed a commercial Lamda40 
spectrophotometer (Perkin Elmer) to acquire sample ex- 
tinction spectra. The sample was an aqueous suspension 
of polystyrene or organosilica microspheres in a square- 
profile (10 mm x 10 mm) quartz cuvette. The concen- 
tration of the suspension was varied (average 1.5 ± 0.6 
g/L) to ensure operation in the single-scattering regime 
without compromising the signal-to-noise ratio. Prior to 
spectroscopy, every suspension was ultrasonicated for 3 
minutes to ensure uniform suspension of microspheres. 
To verify suspension uniformity, spectra measurements 
were repeated immediately after the first recording. Both 
spectra were identical, demonstrating no particle drift 
occurred during measurements. The cuvette was illumi- 
nated by a light source that produced a beam spot of 
approximately 1.0 x 7.5 mm 2 on the front facet of the cu- 
vette, and was characterized by an angular divergence of 
0.4°. A halogen lamp was used as the light source, in the 
visible and near-infrared spectral range, whereas a deu- 
terium lamp was used in the ultraviolet spectral range. 
The illumination light source wavelength was selected 
by means of three spectrally adjacent tunable diffraction 
gratings whose operation wavelength range spanned 300 
- 1100 nm. The detector was a silicon photodetector with 
an area of approximately 49 mm 2 calibrated for the wave- 
length range used. The source to detector configuration 
is shown in Fig. [TJ 

In addition, a cuvette with pure water was placed in 
the reference path of the spectrometer. Optical transmis- 
sion spectra of the sample and reference were recorded 
over the wavelength range of 300 - 1100 nm, and each 
reference spectrum was subtracted from the sample spec- 
trum to account for reflections, aberrations and the in- 
trinsic absorption profile of the solvent and cuvette ma- 
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FIG. 2: Schematic diagram of the optical tweezers refrac- 
tometry for refractive index measurement of individual mi- 
crospheres. A microsphere was trapped by a tightly focused 
laser beam. The microsphere displacement from the centre of 
the trap was detected by illuminating it with a helium-neon 
laser and detecting the signal on the quadrant photodetector. 
The inset shows the microsphere trapped at the centre of the 
focused fiber laser beam. 

The second refractometry method is based on laser 
tweezers 17,[l8|. Here, we measure the particle trapping 
potential created by a high-power focused laser beam, 
as shown in Fig. [21 For a particle displaced by a small 
distance, Al, in transverse direction, the acting force, 
F Q , is measured. From this, the trap stiffness expressed 
in terms of the trap stiffness parameter, a, is calcu- 
lated by F a /Al (measured in piconewtons per nanome- 
tre, pN/nm). This stiffness was calculated for a particle 
of specified size for a range of refractive indices of the par- 
ticle. From the theoretical analysis of measured results, 
as described in Section we obtained the refractive index 
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of the microsphere. In order to create the optical par- 
ticle trap, the beam from a continuous-wave fiber laser 
(power 1 watt, wavelength A = 1070 nm) was focused 
to a diffraction-limited spot using an objective lens with 
high numerical aperture of 1.3. In order to detect the 
particle position in the beam, light from a helium-neon 
laser scattered by the particle was detected on a quadrant 
photodetector. Alternatively, the fiber laser was used 
both for trapping and tracking of the particle position. 
In both cases, the particle displacement from the trap 
centre generated a misbalanced photoelectrical signal on 
the photodetector quantifiable by means of the electronic 
signal processing lj|, 2oj |. 



the ripple structure of the scattering spectrum. 

A local maxima, in Fig. [3)Ja). can be explained by 
the Lorenz-Mie scattering theory. In brief, it is due to 
interference of the incident and diffracted light at the 
particle-medium interface. A phase delay acquired by 
the diffracted light depends on the particle diameter and 
its relative refractive index n sp /n m , where n m is the 
refractive index of the medium. By thin film approxima- 
tion and consideration of the optical path length differ- 
ence, the spectral separation between the maxima of the 
scattering spectrum curve can be determined [ll[ : 



MODELLING METHODS 



AA = (A (m) - A 



A? 



D (n sp - n m ) 



(1) 
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FIG. 3: Typical scattering spectra of microsphere suspensions 
in solution, (a) Fhree distinct features can be observed: rip- 
ples, local maxima, and a sloping curve, (b) Ripple feature of 
the scattering spectra, extracted by high-pass filtering. 

A typical scattering spectrum [Fig. GD^a)] exhibits 
three major features: ripples, local maxima, and a slop- 
ing curve. The ripple represents a small-contrast high- 
frequency sinusoidal signal superimposed on the large 
background signal. An example of the filtered ripple 
structure signal is shown in Fig. E^b). This ripple signal 
is the result of interference between the incident and light 
internally reflected within the particle. It can provide an 
accurate measure of the product of microsphere diame- 
ter D and refractive index n sp . It is important to note 
that the ripple signal is detectable only if the sample is 
monodisperse and the size parameter, x, is much greater 
than 1. Otherwise, the ripple structure of the signal is 
smeared and ultimately undetectable 0, [2l[ . The light 
source and detector configurations have little effect on 



where i and i + 1 stand for the adjacent peaks wave- 
length indices. We note that the fidelity of this rela- 
tionship is affected by the source-detector configuration. 
However, in the context of our work, it provides a good 
estimation of particle refractive index. 

The sloping curve representing the scattering spectrum 
is shown in Fig. [3]Ja) and is tailing towards the long wave- 
length side. This curve also stems from the Lorenz-Mie 
scattering theory. In brief, the optical scatter of a micro- 
sphere takes the form of scattering lobes. The larger 
the size parameter, the more pronounced and narrow 
its forward-scattering lobe. Alternatively, the shorter 
the wavelength of light incident on a sample of micro- 
spheres, the less power is scattered away from the pho- 
todetector, hence the greater intensity is detected. This 
consideration qualitatively explains the sloping trend of 
the recorded spectra, although accurate theoretical mod- 
elling is needed for precise fitting of experimental data to 
the theory, and is beyond the scope of this paper. 

In order to obtain a dispersion curve for the refractive 
index, a Lorenz-Mie scattering model was used to calcu- 
late the extinction cross-section of monodisperse micro- 
spheres suspended in a homogeneous medium. These cal- 
culations assume an incident plane wave and a infinites- 
imally small-size far-field detector i.e. collection angle 
of zero [22j |. Refractive index dependence on the wave- 
length was described by the Cauchy dispersion relation. 
Its general form is given by 



n(X) = A+^ + ^ (2) 

where the Cauchy coefficients A, B (/jm 2 ) and C (/im 4 ) 
are unique to the material [23j |. Using an unconstrained 
nonlinear optimization method implemented in Matlab, 
particle size and a dispersion curve of refractive indices 
over a wavelength range of 300-1100 nm were deter- 
mined. 



Optical tweezers refractometry 

Optical tweezers refer to a technique of trapping a sin- 
gle microsphere immersed in fluid at the focus of a laser 
beam. It is well known that particles in fluid experi- 
ence random Brownian motion. For a trapped micro- 
sphere, Brownian motion is suppressed, where the sup- 
pression efficiency is a monotonic function of the optical 
trap stiffness, a, which in turn depends on the micro- 
sphere size and refractive index. Therefore, the analysis 
of the Brownian motion amplitude versus the controllable 
stiffness provides a basis for the determination of the par- 
ticle refractive index. The Brownian motion amplitude is 
determined by detecting the microsphere position fluctu- 
ation in the trap, which is analysed in terms of its power 
spectrum amplitude. The roll-off frequency, / , of this 
spectrum provides a reliable measure of the trap stiffness 
via the following relationship 17 1: 




FIG. 4: SEM image of monodisperse organosilica micro- 
spheres 



fo = 



2irf3 



(3) 



where (3 is the drag coefficient well tabulated for var- 
ious fluids, a relates to x and n sp . This dependence 
has been analysed in the framework of the generalized 
Lorenz-Mie theory and detailed elsewhere [l7], 24, (2||. 



Experimental measurement of a and its comparison with 
the theoretical curve of a versus n sp at a fixed x, yields 
the refractive index of the particle. 



RESULTS AND DISCUSSION 

A low-dispersion particle sample is required for pre- 
cise determination of refractive indices of microspheres by 
means of scattering spectroscopy. Fig. [4] shows a typical 
ensemble of organosilica microspheres imaged by SEM. 
The CV of organosilica microspheres was determined to 
be less than 10%, suitable for the scattering spectroscopy 
measurements. 

Organosilica microspheres with well controlled diam- 
eters could be produced in the size range of 1-6 /im. 
In Fig. G^a), (b), and (c), recorded and calculated 
spectra for 2.01 and 5.26 fim polystyrene and 5.94 /xm 
organosilica microsphere samples, respectively, are pre- 
sented. Large size organosilica microspheres were cho- 
sen to ensure the experimental spectra exhibited at least 
two maxima. In the case of the monodisperse 2.01 /zm 
polystyrene microsphere aqueous suspension sample, the 
recorded scattering spectrum was accurately fitted by the 
calculations using the Lorenz-Mie theory. The calcula- 
tions of the dispersion of polystyrene refractive index are 
presented in Fig. [UJa). The theoretical fitting curve was 
parameterized in terms of the Cauchy coefficients (pre- 
sented in Table HI, which shows very good agreement with 
those reported in the literature. We believe that the good 



agreement between the experimental and theoretical data 
was due to the relatively small sphere size, which forward- 
scattered incident light over wide angles. The evaluation 
of the 2.01 /zm polystyrene sphere scattering showed that 
the forward-scattering lobe overfilled the area of the pho- 
todetector for A < 390 nm, as illustrated in the Fig. [SJa) 
inset. This qualitative model also explains the discrep- 
ancy between theoretical and experimental data at A < 
390 nm, where the next side-lobe incident on the detector 
interferes constructively with the original incident plane 
wave, unlike the forward-scattering lobe which interferes 
destructively (which is why extinction in the forward di- 
rection occurs). Therefore, the extinction is reduced in 
comparison with that calculated in the small photode- 
tector approximation. Alternatively, one can view this 
in terms of the size parameter x. In Fig[5]Ja), the intro- 
duction of the second side-lobe can be observed at an x 
value of 17. The effect of multiple side-lobe interference 
is much more pronounced in the case of the scattering 
spectrum of 5.26 /im polystyrene microspheres, shown in 
Fig. EJb), and 5.94 /im organosilica microspheres, shown 
in Fig. [5Jc) . Here, the maximum x values are 15 and 
17. The insets of Fig. UK a) and (b) illustrate the detec- 
tion scenarios for 2.01 and 5.26 /im polystyrene sphere 
aqueous suspension, respectively. Obviously, reducing 
the detector size can avoid the collection of the side- 
scattering lobe. However, commercially available spec- 
trometers have relatively large size detectors, which are 
designed to maximize detection sensitivity. 

For an accurate theoretical model, which will result in 
precise fitting of measured data, a finite detector must 
be accounted for. However, this will require more com- 
plex Mie theory calculations and enormous computing 
power. Therefore to obtain maintain the simplicity of 
this method, we employ high pass filtering to isolate the 
ripple signal from the sloping curve of the recorded spec- 
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FIG. 5: Recorded and calculated spectra of microsphere suspension of (a) 2.01 fj,m and (b) 5.26 /im, polystyrene and (c) 5.94 
/xm organosilica microspheres in water. High-pass filtered spectra of (a') 2.01 /im and (b') 5.26 /im polystyrene, and (c') 5.94 
/jm organosilica microsphere suspensions. 
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FIG. 6: Refractive index dispersion of (a) polystyrene and 
(b) organosilica microspheres by scattering spectroscopy (SS) 
and laser tweezers refractometry (LTR). "Ma et al." denotes 
for polystyrene, as given by Ma et al. (2003) p]§]; "Malitson" 
denotes refractive index dispersion for fused silica, as given 
by Malitson (1965) [H. 



tion, was not accounted for in our modelling method. 
The effect of multiple scattering has been shown, in liter- 
ature, to increase with optical depth (concentration) and 
gradually smoothes out the maxima and minima [27] |. 
The spectra of solutions with increasing concentrations 
were recorded, and did not exhibit the effects of multi- 
ple scattering, which were observed in the work of Cohen 
(1975). 

The refractive index of the microspheres was also 
measured using optical tweezers refractometry. This 
method did not require monodisperse suspensions since 
only a single particle was measured at a time. For 2.09 
firm polystyrene and 5.16 /jm organosilica microspheres, 
the power spectrum amplitude of Brownian motion, 
shown in Fig. , was recorded, and the roll-off frequencies 
determined. Using Eq. ^ and tabulated drag constants, 
the trap stiffness was calculated, and is shown in Table 
HIl By using experimentally measured diameters of 
the trapped polystyrene and organosilica microspheres 
we calculated a versus n sp in the framework of the 
Lorenz-Mie scattering theory, with the results presented 
in Fig. [8] The intersection of the theoretical curves and 
corresponding horizontal lines of the determined values 
of a yielded the refractive index of the trapped mi- 
crospheres. For polystyrene microspheres, we obtained 
a refractive index of n sp = 1.58 ± 0.01 at A — 1070 
nm. This value deviates by a 0.6% from that given in 
literature [l3j]. For the organosilica microspheres, we 
obtained refractive indices of 1.52 ± 0.01 and 1.51 ± 
0.01 at 589 and 1070 nm respectively, which was in a 
good agreement with the respective values of 1.518 and 
1.511 obtained by the scattering spectroscopy method. 
For comparison with the scattering spectroscopy results, 
the refractive index values obtained by laser tweezers 
refractometry are also shown in Fig. [BJa) and (b). 



tra. The results are presented in Fig. EJa') and (b') for 
2.01 and 5.26 /im polystyrene microspheres, respectively. 
For the 5.26 /im polystyrene sample, we used the ripple 
signal to determine the refractive index dispersion rela- 
tionship, provided an independent measurement of the 
particle size could be made. For the organosilica micro- 
spheres, the ripple structure was not observable due to 
relatively large CV value. Therefore, we used the local 
maxima to determine their refractive index. The max- 
imum was isolated by a high-pass filter to suppress the 
sloping feature of the recorded spectrum. The resultant 
spectrum is shown in Fig. Etc') , and displays a reason- 
able agreement with the theoretical scatter spectrum. 

Multiple scattering, other than in the forward direc- 



CONCLUSIONS 

The refractive index of organosilica microspheres was 
determined in the 300-1100 nm wavelength range using 
laser tweezers refractometry and a relatively simple scat- 
tering spectroscopy technique. The refractive index of 
organosilica has been found to range from 1.60-1.51 over 
the wavelength range of 300-1100 nm, which is consider- 
ably greater than that of the organosilica substrate, fused 
silica. The development of the refractometry methods 
addressed in this paper can be useful for accurate deter- 
mination of nano(micro)material refractive index in situ, 
whereas the tabulation of refractive index of organosil- 
ica microspheres has important applications to biological 
screening using flow cytometry. 
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TABLE I: Cauchy coefficients for refractive index dispersion relation (Eq. [2} 
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FIG. 7: Plot of power spectral amplitude densities of 2.09 
fim polystyrene (PS) and 5.16 ^tm organosilica microspheres 
(ogSI) versus frequency of particle Brownian motion in the 
optical tweezers trap. The roll-off frequency (f ) was used to 
determine the trap stiffness. 



TABLE II: Parameters determined by Laser Tweezers Refrac- 
tometry 
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